Quantum computing with an inhomogeneously broadened ensemble of ions: 
Suppression of errors from detuning variations by specially adapted pulses and 

coherent population trapping. 
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The proposal for quantum computing with rare-earth-ion qubits in inorganic crystals makes use 
of the inhomogeneous broadening of optical transitions in the ions to associate individual qubits 
with ions responding to radiation in selected frequency channels. We show that a class of Gaussian 
composite pulses and complex sech pulses provide accurate qubit 7r-rotations, which are at the same 
time channel selective on a 5 MHz frequency scale and tolerant to ±0.5 MHz deviations of the 
transition frequency of ions within a single channel. Rotations in qubit space of arbitrary angles 
and phases are produced by sequences of 7r-pulses between the excited state of the ions and coherent 
superpositions of the qubit states. 

PACS numbers: 03.67.-a, 33.25.+k, 82.56.Jn 



I. INTRODUCTION 



Rare-earth ions doped into inorganic crystals experi- 
ence inhomogeneous broadening. If the doping concen- 
tration is high, a near field probe is able to address in- 
dividual ions and to induce resonant excitation hopping 
among nearby ions which can hereby be used for quan- 
tum computing lj. Along similar lines, the rare-earth- 
ion-doped inorganic crystal proposal for quantum com- 
puting (REQC) was developed Instead of excitation 
hopping, it uses a dipole blockade mechanism, similar to 
the Rydberg blockade mechanism for neutral atoms ||, 
which conditions gate operations on the interaction be- 
tween a pair of ions exceeding a certain threshold rather 
than having a precise value. The REQC proposal iden- 
tifies qubit registers as entire ensembles of collections of 
interacting ions, selected by their response to radiation 
at certain frequencies, and distributed over the spatial 
extent of the crystal. 



There is a trade-off between the number of registers 
accepted by the spectroscopic selection and the preci- 
sion we can request for their transition frequencies and 
coupling strengths, and hence, ultimately, between the 
scalability of the concept and the fidelity of the gate op- 
erations. Gate operations are needed that can be carried 
out with high fidelity, even if the physical parameters 
of the system are allowed to vary, and following strong 
and convincing recommendations by Jones jj], we apply 
composite rotation techniques from liquid and solid state 
NMR to the problem. Our method, which also involves 
quantum optical interference techniques, is described in 
detail for the REQC scheme, but in the conclusion we 
shall comment on its possible application to a number of 
other quantum computing proposals. 



II. RARE-EARTH-ION QUANTUM 
COMPUTING SCHEME 



In the rare-earth ion scheme for quantum computing, 
qubits are encoded in two ground hypcrfmc levels |0)and 
|1) of the dopant ions, which are only negligibly disturbed 
by the interaction with the crystal host. Excited states of 
the ions interact more strongly with the crystal, and this 
gives rise to significant broadening of the optical transi- 
tion between the qubit states and the excited state. The 
characteristic feature of the proposal is the identifica- 
tion and initialization of the qubits and qubit registers 
in the crystal: The ith qubit is chosen as the ensemble of 
ions within the inhomogeneously broadened line which 
have their absorption frequency at a chosen frequency 
Vi . Ions in a frequency window around the qubit channel 
are pumped to auxiliary states to prevent them from be- 
ing excited by the near resonant laser light at the qubit 
frequencies, see Fig. ^ 

The excited state of an ion establishes a surrounding 
permanent dipole field, which shifts the transition fre- 
quencies in neighboring ions and which can hence con- 
trol their interaction with a laser field. Pairs of interact- 
ing qubits are hence identified, or rather selected, 
spectroscopically, by pumping away those ions which do 
not have their transition frequency shifted out of the fre- 
quency window around the channel frequency VjU) by the 
dipole-field, when neighboring ions in the channel are 
excited. Physically, the crystal still contains all the ions, 
but around frequency channels VjU), the ions, populating 
the qubit levels, are found in close lying pairs, since all 
ions with insufficient excited state interaction have been 
deported to passive spectator levels. The ensemble of 
ions in frequency channels Vi and Uj thus constitute pairs 
of mutually interacting ions. Theoretical estimates Q 
and preliminary experiments 5] suggest that the dipole 
energy shifts easily amount to several MHz, and the ions 
can be excited with Rabi frequencies in the MHz range. 
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FIG. 1: A schematic diagram of a part of the inhomogeneously 
broadened spectrum of excitation frequencies to the excited 
state |e) in rare-earth-ions in an inorganic crystal. Two qubits 
i and j are defined by two different transition frequencies in 
the spectrally hole-burnt structure. 
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FIG. 2: A C-NOT operation with qubit i as control and 
qubit j as target bit. The operation makes use of the fol- 
lowing sequence of 7r-pulses: 1. |0} i *-* \e) i , 2 



explanation. 



10), 



|e>., and 5. |0>. «-> \e) i 



|0>,-|e)„3. 
See text for 



that the transitions 1-5 in Fig. will not be resonant 
for all ions, and this seriously compromises the gate fi- 
delity in the system. One could apply very short pulses 
which are less sensitive to the resonance criterion, but 
they would excite ions outside the windows around the 
frequency channels, and they would not display the re- 
quired sensitivity to the frequency shift on ion j due to 
the excitation of ion i in the two-bit gate. In a recent 
paper Q it was proposed to use a combination of adia- 
batic processes and coherent control theory to increase 
the tolerance to variations in physical parameters in a 
quantum computing model, but very high fidelity was 
only achieved with very slow processes, which would be 
compromised by decoherence and decay in the REQC 
proposal. Instead, we shall present and compare the ad- 
vantages of two different realizations (composite rotation 
and complex hyperbolic secant pulse) of 7r-pulses, which 
rapidly transfer population completely from one state to 
another in a two-level system, as long as the detuning 
of the transition lies within a 1 MHz wide interval, but 
which does not transfer any population if the detuning 
exceeds just a few MHz. Composite pulses were, in fact, 
already used in quantum computing experiments with 
trapped ions [jj, not to correct for variations in phys- 
ical parameters, but to obtain it and 4tt Bloch sphere 
rotations on two different transitions where the Rabi fre- 
quencies differed by a factor of We shall propose 
a combination of coherent population trapping and 7r- 
pulse rotations as a means to make robust qubit rotations 
through arbitrary angles and on arbitrary qubit states. 



A. Composite pulses 



The homogeneous linewidth is small corresponding to ex- 
cited state lifetimes on the msec scale, whereas the qubit 
states have lifetimes of seconds with coherence times in 
the msec range with possibilities for increases up to sev- 
eral hundred msecs in the presence of magnetic fields [(| . 

The implementation proposed in Ref. ,2] for a C-NOT 
operation, with qubit i as control bit and qubit j as target 
bit, is shown in Fig.|3 If qubit i is initially in state |0)i, 
it is excited by the first pulse, and the resulting frequency 
shift of qubit j makes the steps 2, 3 and 4 non-resonant, 
hence nothing happens, and qubit i is finally returned 
to its initial state. If qubit i is in state there is no 
frequency shift, and the resonant processes 2, 3 and 4, 
effectively exchange states |0).,- and as desired for 
the operation. 



III. HIGH-FIDELITY tt-PULSES 

To have a sufficiently large number of active ions, and 
hence an appreciable number of quantum registers with, 
e.g., three or more coupled qubits, it is necessary to use 
ions with frequencies in a not too narrowly defined inter- 
val. The widths of the qubit channels, however, imply 



Composite pulses can be assigned to two main classes, 
A and B, where pulses of type A produce a fully compen- 
sated rotation of the system for all initial states, whereas 
pulses of type B provide the compensated transformation 
only for particular initial states. Levitt has tabulated 
several sequences of rectangular pulses that provide net 
90° and 180° Bloch sphere rotations of two- level systems 
and which are more tolerant to frequency detunings than 
a single pulse. Because of the high frequency components 
of the rectangular pulses, these pulses do not, however, 
satisfactorily exclude excitation of the surrounding ions, 
and we have instead studied composite pulses based on 
Gaussian field envelopes. On resonance, the accomplish- 
ments of a laser pulse is governed only by pulse area and 
phase, and hence it is natural to replace a rectangular 
pulse sequence with a sequence of Gaussian pulses with 
the same areas and phases i-e., with the time 

dependent Rabi frequencies 

n R (t) = V e j^l e -{t-t,)V^\ t G [tk Mfc + a] . 

(1) 

Composite rotations correct errors to high order, and 
there is a good portion of experience and magic in the 
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FIG. 3: Excitation by a composite pulse sequence 
909ol80i8u909o of rectangular pulses (left) and a compos- 
ite pulse sequence 92.50g6.98 192. 006. 8692. 42g6. 23 of Gaussian 
pulses (right). The total duration of both sequences is 1.5 /xs. 



design of good sequences. Merely replacing rectangular 
pulses with Gaussian pulses of equal pulse area and phase 
does therefore not optimally preserve the error compen- 
sating properties. The composite rotation with rectan- 
gular pulses, 90gol80 90 9 o [10|. provides an error com- 
pensated 7r-pulse of type B, and the left panel in Fig. [3| 
illustrates its tolerance to detuning errors on the 1MHz 
scale, but also its erroneous excitation of ions more than 
5MHz away from the channel frequency. Gaussian pulses 
with the same areas give a better performance without 
the off-resonant excitations, but we may do even bet- 
ter by slightly adjusting the pulse areas and phases. We 
have done this in a variational approach, assuming three 
pulses, and hence six free area and phase parameters. For 
comparison, we choose the same duration of the rectan- 
gular pulse and the Gaussian pulse sequences and a tem- 
poral cut-off a = 3.5cr, so that the pulses are truncated 
when the amplitude is 0.2% of its maximum. By insisting 
on perfect operation on resonance and minimum errors 
for a set of detunings, we have used a variational search to 
identify an optimum Gaussian pulse sequence with pulse 
areas and phases 92.50g6.9sl92.006.8692.42g6.23. Simu- 
lations of the excitation with this composite pulse are 
plotted in the right panel of Fig. |3J It has been been 
shown [llj that composite pulses can satisfactorily deal 
with errors in coupling strengths, e.g., due to spatial in- 
homogeneities of the field strength or variations of the 
transition dipole moments in the ions, but a gate that 
works for all input states was not identified in that work, 
and no solution was presented for the important case of 
detuning errors in the REQC proposal. 



B. Complex hyperbolic secant pulse 

As shown by Silver et al. |l2j |. complex hyperbolic se- 
cant pulses offer an alternative process, of class B, with 
the desired properties. The complex Rabi frequency, 



FIG. 4: Excitation by a complex hyperbolic secant pulse for 
different detunings. The left panel shows the time evolution 
on the Bloch sphere of ions with different detunings. The 
right panel shows the excited state population after the pulse 
as function of detuning. The parameters are n = 3, fin = 2 
MHz and f3 — 0.64 MHz, and the duration of the pulse is 3 fis. 



of the complex hyperbolic secant pulse is given 

by 



n R (t)=Cl {sech[/3 (t - t )}} 



l+i/i 



(2) 



where [i is a real constant, Qq is the maximum Rabi fre- 
quency and P is related to the FWHM of the pulse. This 
corresponds to a real sech pulse envelope in combination 
with a tanh frequency sweep around the carrier laser fre- 
quency vl'- 



\n R {t)\ = n sech (13 (t~ t Q )) 
v(t) = i/£ + /i/3tanh(/3(t- 



to))- 



(3) 
(4) 



The excitation by such a pulse is shown in Fig. 0] as a 
function of the detuning (v^ — v£)- In addition to its 
favorable tolerance to detuning errors, the complex sech 
pulse is independent of pulse amplitude if only it exceeds 
a critical threshold. Its disadvantage compared to the 
composite Gaussian pulse is the double implementation 
time. 



IV. ARBITRARY ROTATIONS, 
COMPENSATION OF PHASE SHIFTS 

A. Using dark states for arbitrary qubit rotations 

The pulses described are excellent for steps 1-2 and 4-5 
in the C-NOT scheme in Fig. |2 since the initial states 
are explicitly known for all these processes due to the 
fact that the excited state is unpopulated when the op- 
eration starts. The initial state before step 3, however, 
is an unknown superposition state of and |e)j, and 
a class A pulse is required to perform a perfect 7r-pulse 
on that state. Such composite pulses exist based upon 
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square pulses, but they are difficult to build with Gaus- 
sian pulses, and, we also find it difficult to make compos- 
ite pulses for arbitrary rotations, since the ones provided 
in the literature often involve a larger number of pulses of 
which some are rather short, and hence, they will interact 
with ions outside the hole burnt structures. Instead, we 
propose an alternative technique, which does not require 
class A pulses but nevertheless offers robust arbitrary ro- 
tations in the qubit space. To this purpose, we shall make 
use of the full three-level structure of the ions and turn 
on fields with complex Rabi frequencies £ljio(t)e~ lVo and 
Qm{t)e~ lipi that simultaneously couple \0)j and to 
\e)j. The system now has two ground state superposi- 
tions, that will be respectively coupled and uncoupled to 
the excited state. If Clno(t) = ^m(t) we can define a 
new orthonormal basis according to 



(|0> 



1 

V2 



(5) 



If the relative phase of the two fields are chosen so that 
fi ~ fo = 0j only |1) couples to the excited state while 
|0) is a dark state, unchanged by the interaction with 
the fields, and vice versa if ipi — tpo = <ft + ir. Apply- 
ing constant phase factors e~ llpo , e~ lipi on the sequence 
of Gaussian pulse or the complex sech pulse described 
above, we can therefore implement a robust 7r-pulse of 
class B, i.e., from the south pole to the north pole on the 
Bloch sphere of states |e) and |T). If we drive the state 
back down again by a new 7r-pulse with a phase that is 
shifted by 7r + 9 in comparison to the first one, the net 
effect becomes a robust phase shift on state |T), 
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(6) 



In the original qubit basis (|0),|1)), this procedure is 
equivalent to the unitary operation 



U = e* 



cos 



2 

sin ■ 



ie** sin ! 
cos£ 



(7) 



i.e., a general unitary rotation in the qubit state basis. A 
NOT operation is achieved if 9 and (j> are both set to ir, 
but we note that this approach offers any arbitrary qubit 
rotation, based exclusively on class B 7r-pulses which, in 
turn, can be implemented in a robust manner by the 
pulses analyzed in Figs. and A similar combination 
of coupled and uncoupled superposition states and the 
STIRAP process was recently proposed [l3|. 
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TABLE I: Implementation scheme for a robust C-NOT op- 
eration. Pulses 4 and 5 compensate unknown phase factors, 
dependent on the detuning of the target bit, caused by pulses 
2 and 3. Pulses 7-12 compensate phase factors dependent 
on the detuning of the control bit, caused by the sequence 
1-6. All 7r-pulses should be implemented either as composite 
pulse sequences or complex sech pulses to yield the desired 
tolerance to detuning variations. 



same phase to be acquired by |0). We hence suggest to 
apply compensating pulses, which have no other effect 
than to introduce the same detuning dependent phase 
shift on that state, that is a composite or complex sech 
7r-pulse followed by the inverse 7r-pulse, i.e., with a phase 
shift of 7r (but no 9), on |0) <-> |e). The choice of im- 
plementation and the duration of these pulses must be 
the same as those applied to |1) <-> |e) for the detuning 
dependent phases to be accurately compensated. 

If the operation is a C-NOT operation, i.e. a two-bit 
operation, we will come across two different detuning de- 
pendent phase factors, one dependent on the detuning of 
the target ion and the other dependent on the detuning of 
the control ion. Since the detunings, and hence the phase 
factors, are unknown and not necessarily equal, we have 
to make sure that every qubit state acquires both phase 
factors in order to be able to consider the total unknown 
phase as global. The target qubit rotations, described 
above, only took place for the state vector components 
with the control qubit in the state |1), and, thus, condi- 
tioned on the control bit being in state |0), we apply a 
new sequence of compensating pulses to the target states 
|0) and |D so that all states |00), |01), 1 10) and |11) finally 
are equipped with the same detuning dependent phases. 



B. Phase compensation 

The Bloch sphere rotations introduce detuning and 
time dependent phases on the individual ions. A phase 
factor can be considered as global, and thus disregarded 
only if it appears in front of all qubit states. The two 
states involved in the above rotations are |1) and |e) of 
the target bit, and one must necessarily arrange for the 



V. COMPLETE GATE SEQUENCE, 
NUMERICAL RESULTS. 

The implementation scheme for a C-NOT operation 
with qubit i as control bit and qubit j as target bit is sum- 
marized in Tabled and we draw the attention to the fact 
that 12 7r-pulses (which are all composite or complex sech 
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FIG. 5: Simulations of a C-NOT operation performed with 
optimized composite Gaussian pulses (left) and complex sech 
pulses (right). The initial qubit state was ip = \/0.1 1 00) + 
\/fl2|01) + \/O3|10) + V0A\11). The dotted vertical lines 
represent the edges of the qubit channels and of the spectral 
hole, respectively. The total duration of the operation is 18 fis 
if all 7r-pulses are implemented with composite pulses and 
36 /xs if they are implemented with complex sech pulses. 



pulses, but not of class ^4) are applied. This scheme also 
suffices if the desired gate is a given rotation of the target 
bit, conditioned on the control bit; only the phase of pulse 
3 in the sequence should be adjusted accordingly. Simu- 
lations according to this scheme have been performed by 
solving the equations of motion in the space spanned by 
the nine states \ciibj), where a,b — 0,1, e. An example 
of the final state after action of the C-NOT operation 
on an initially entangled state of the qubits is shown in 
Fig. For simplicity, the detuning parameter was as- 
sumed to be the same for the control and target ions in 
the calculations leading to Fig. [S] The C-NOT operation 
effectively exchanges the amplitudes and hence the pop- 
ulations on states |10)and |11), and this is precisely what 
we observe for all detunings in the central 1 MHz wide 
frequency interval, whereas nothing happens to the state 
of ions more than 5 MHz away from the exciting lasers. 
The intermediate frequency intervals are irrelevant, since 
they contain no ions, see Fig. ^ Decay and decoherence 
become important on the msec time scale, and it is hence 
important that the gates are fast on this time scale. The 
total duration of the operation is 18 /is if all 7r-pulses are 
implemented with composite pulses and 36 /is if they are 
implemented with complex sech pulses. 

In Fig. H3 we show the results of our simulations when 
the control and target detunings are allowed to both vary 
independently within the central 1 MHz wide frequency 
channel. The four sub-plots show the populations, which 
are also presented in Fig. [3 and we see that the devi- 
ations on the edges of the frequency interval are barely 
10~ 3 when the 12 7r-pulses are realized with complex sech 
pulses whereas implementation with composite Gaussian 
pulses results in deviations on the order of 10 -2 . The 
relative phases between the four states are accurate to 
the level of 1 degree for both kinds of pulses. 
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FIG. 6: Simulations of a C-NOT operation where the con- 
trol and target detunings are allowed to vary independently 
within the qubit frequency channel. The initial state was 
ip = VOT i00) + ^O2!01) + V03 1 10) + V04 and (a) is 
realized with composite Gaussian pulses while complex sech 
pulses are used in (b). Note the different scales. 



VI. DISCUSSION 

We have shown that composite and complex sech 
pulses can be applied in combination with coherent pop- 
ulation trapping to perform accurate quantum gate op- 
erations on ions in the REQC quantum computing pro- 
posal. The composite Gaussian pulses do not reach the 
same high fidelity as the complex sech pulse. But, they 
are twice as fast and higher accuracies may be achieved 
with improved optimization. The combined tolerance to 
detuning errors on one frequency scale and frequency se- 
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lectivity on a single order of magnitude larger scale are 
necessary to make quantum computing feasible at all in 
the REQC proposal, and the first experiments to test 
our gates are in preparation. Our theory also provide 
important ingredients for the further development of the 
proposal, in particular it contributes significantly to the 
scalability to larger numbers of coupled qubits. 

The achievements of the fault tolerant pulses are re- 
markable, and it is clear that many quantum optics and 
quantum information tasks may benefit from their intro- 
duction. Other proposals exist for quantum computing 
where individual addressing of qubits is made spectro- 
scopically, e.g., of optically trapped neutral atoms in in- 
homogeneous magnetic fields , of trapped ions in in- 
homogeneous magnetic 0] and optical dipole fields . 
Here, tolerance to variations in the detuning may be de- 
sirable as in the present work. Likewise, systems which 



use spatial addressing may suffer from difficulties of fo- 
cusing the fields on single qubits, and here tolerance of 
the relevant qubit to variations around the maximum 
field and of the neighboring qubits to small variations 
around the desired vanishing field can be accommodated 
by only slight modifications of the techniques presented 
in this paper. 
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